Cavity-QED tests of representations of canonical commutation relations employed in 

field quantization 



> 



23 



Marek Czachor and Marcin Wilczewski 
Katedra Fizyki Teoretycznej i Metod Matematycznych 
Politechnika Gdanska, Narutowicza 11/12, 80-952 Gdansk, Poland 

Various aspects of dissipative and nondissipative decoherence of Rabi oscillations are discussed in 
the context of field quantization in alternative representations of CCR. Theory is confronted with 
experiment, and a possibility of more conclusive tests is analyzed. 

PACS numbers: 42.50.Pq, 42.50.Xa, 03.70.+k 



I. CAVITY QED IN DIFFERENT REPRESENTATIONS OF CCR 



o ■ 

o : 

Although the notion of entanglement between atomic and electromagnetic degrees of freedom plays a central role 
in quantum computing architecture based on cavity QED [lll^lSl^|a,lajtD)l3jl2() the very concept of entanglement 
leads to conceptual difficulties if quantum vacuum comes into play (cf. uniqueness of the vacuum versus violation of 
the Bell inequality |Toj . problems with teleportation of quantum fields ambiguous entanglement with vacuum 
CN| \ [l2l ^3)- One of the problems is that the electromagnetic field can be quantized in different representations of 
■ canonical commutation relations (CCR). As shown in [jjj the degree of entanglement is a representation-dependent 
property, and it is not clear which representations are really physical. The problem is a part of a wider and ongoing 
\l discussion on different quantization paradigms [l4|. 

Now, can the available experimental data distinguish between different representations of CCR? The answer is less 
obvious than one might expect. In this Letter we will try to clarify the status of some data from cavity QED, and 
y—{ ■ then discuss possibilities of more definitive tests. 

| We first analyze at a representation independent level the simple problem of Rabi oscillation of a two-level atom in 
an ideal cavity (for technicalities we refer to |l5|). In the second step we take into account two types of decoherence 
that should occur in realistic experiments. Following [Tfil ] we distinguish between dissipative and nondissipative 
decoherence and model dissipation employing the results of Then we compare theoretical predictions based on 
irreducible representations with the experimental data of the Paris group |l8|. Our conclusions are basically consistent 
Qh! with both and 0]: The observed decoherence appears to be entirely of a nondissipative type, but it is not clear 
(Jj ■ why the effect of dissipation is invisible. Perhaps the fact that a photon is with probability 1 absorbed by the atom 
at times separated by the Rabi period leads to a sort of Zeno effect. This point requires further experimental and 
theoretical studies, and is beyond the scope of the present paper. 

Assuming that Rabi oscillations indeed do not reveal observable damping due to energy dissipation we ask to what 
extent the experiment can distinguish between reducible and irreducible representations of CCR. In physical terms 
the question can be translated as follows: How many oscillators do we need to model quantum fields? The standard 
answer is that we need one oscillator per mode. We show that in reducible representations the data only set certain 
limitations on the number of oscillators, and this number is independent of the number of modes. 

Finally, we suggest that one should repeat the measurements reported in |18| with better cavities, finer time 
resolution, and monitor the Rabi oscillation for longer times. The point is that the decay due to experimental 
imprecisions may mask quantum beats of a completely new type and origin. In principle, the beats can be observed 
in a form of vacuum collapses and revivals, the effect occurring in reducible A-representations |19|. Observation of 
the revival would be of fundamental importance for our understanding of field quantization. 



II. RABI OSCILLATIONS AT A REPRESENTATION-INDEPENDENT LEVEL 

Similarly to (lsf we work with the Jaynes-Cummings model poll2l|. The crucial point is that we begin with solving 
Heisenberg equations of motion for the two-level atom at a representation independent level j^. The CCR algebra, 
in its general form, reads 

[a k ,a* k ,] = Skk'h, (1) 

Ik commute with all the other operators, and I£ = Ik- We do not assume that Ik is proportional to the identity 
(this generality will pay, as we shall see shortly). By Schur's lemma Ik is necessarily proportional to the identity 
only in irreducible representations. We employ the usual notation 21] where Ri = er;/2, R± = R\ ± 1R2, ai are the 
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Pauli matrices, and g is a complex coupling parameter. We assume there exists a free-field Hamiltonian Hq satisfying 
[ak,H Q ] — u> k a k , [ a *kiH Q ] = —u) k a* k . Note that H cannot, in general, be given by ~}2, k Uka k a k ; the latter works only 
for some representations (e.g. for irreducible representations with I k equal to an identity, or for the reducible 'TV = V 
representation; 'TV > V reducible representations require a different construction). 

Let us now select a frequency tu p = \p\ = u> and assume that only this frequency couples to the two-level system. 

We also split Hq into two parts: H^ commuting with a p and a*, and Hq — u>N p , where [a p , N p ] = a p , [a*, N p ] = —a*. 
The model is given by the full Hamiltonian 

H = ui R 3 + H a + gR+a p + gR~a* p . (2) 

Solving the Heisenberg picture equations we find 



a sm^OfltU fA sm 2 (Q R t) . sm(2Q R t) \ fA sm\Q R t) . sm(2Q R t) ■ „ 



(3) 



where A = w — f2,R = \/A 2 /4: + \g\ 2 X, and 

X = {R 3 + l/2)I p + a* p a p . (4) 
The next important notion that can be introduced at a general level is the displacement operator 

D(z) = exp^ (z k a* k - z k a k ). (5) 
fc 

Acting with D(z) on a vacuum vector we obtain a coherent state. Its form depends on what is meant by vacuum in 
a given representation. 

We will not discuss in more detail the irreducible representations since, as shown in |lfij |. they all yield physically 
equivalent and well known results. Instead, we directly turn to the 'N < oo' reducible representation introduced in 
|23| and worked out in many details in [2J, [25|, |2(j . 



III. TV < oo REPRESENTATION 



The representation is constructed as follows. For simplicity we ignore here the polarization degree of freedom (see 
however |2J1 l24l 123. |2jj). Take an operator a satisfying [a, a*] = 1 and the kets |fc) corresponding to standing waves 
in some cavity. We define 

a k = \k){k\<8)a, h = \k){k\ ® 1. (6) 

The operators © satisfy (JTJ, where 8 kk i is the 3D Kronecker delta. The fact that I k is not proportional to the 
identity means that the representation is reducible. In our terminology this is the 'TV = 1 representation'. Its Hilbert 
space TL is spanned by the kets \k,n) = \k)\n), where a*a\n) = n\n). Such a Hilbert space represents essentially 
a single harmonic oscillator of indefinite frequency (for physical motivation cf. [23L |24| and the Appendix in [l5jp. 
An important property of the representation is that Ik = I is the identity operator in H. A vacuum of this 
representation is given by any state annihilated by all a k . The vacuum state is not unique and belongs to the 
subspace spanned by |fe,0). In our notation a TV = 1 vacuum state reads \0) = ^2 k O k \k,0) and is normalized by 
y\ \O k \ 2 = Y^k ^ k ~ ■*-> Z k = |Ofc| 2 . Such a vacuum represents a single-oscillator ground-state wavepacket. As shown 
in |2 ll l25j in a fully relativistic formulation the maximal probability Z = max k {Z k } is a Poincare invariant and plays 
a role of renormalization constant. For N > 1 the representation space is given by the tensor power 7t = Tt® N , i.e. 
we take the Hilbert space of TV (bosonic) harmonic oscillators. Let A : TL — > TL be any operator for N = 1. We denote 
A (n) = j®(n-i) A 7 ®(JV-n) ; A (n) . 21 for 1 < n < TV. For arbitrary N the representation is defined by 

N N 
n—l n—1 

im,,aZ>] = s kk 'h, z2h=L = i® N (8) 
k 
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and the iV-oscillator vacuum is the TV-fold tensor power of the TV — 1 case, a kind of Bose-Einstein condensate 
consisting of TV wavepackets: 

\0) = \0)®---®\0) = \0)® N . (9) 
The free-field Hamiltonian is, for TV = 1 and uj k = |fe|, 

H a = Y2 UJ ka* k a k = ^^uj k \k){k\® a*a. (10) 

fc fc 

In each eigensubspace with fixed |fc) the operator Hq is just an ordinary Hamiltonian of the oscillator with frequency 
Lu k . Let us note that one can work also with 

i ^ 

H = ^^2uJk{a* k a k + a k a* k ) = qjfcjfc) <fcj <g) a* a + - ^uj k I k . (11) 

fc fc fc 

The vacuum term is a well-defined Hermitian operator, and can be removed by a well-defined unitary transformation. 
This is an example of a procedure that can be termed, after Finkelstcin, "regularization by quantization" 14] . For 
arbitrary TV the generator of free field evolution is the Hamiltonian of TV noninteracting oscillators, i.e. H_q = 

X!^Li Hq U \ Let us stress that H_ should not be confused with k-'fcOfcOfc- The operator a k a k nevertheless occurs in 
(|4*|) and thus plays an important role in the Jaynes-Cummings problem. Our definition of H_ implies that [a k ,H_ ] — 
ui k a k which is the formula we required at the representation independent level. 
A monochromatic coherent state with frequency lo is given by the usual formula 

\z) = exp (za* - za p )\0). (12) 

Starting with the excited state and a vacuum field, \^f) — \+}\0), we find 



w(t) - WMW) ~ 2 " 2^ \9\ A? Ay4+\ g \z S /N { a ) ^ (1 " Zp) ' (13) 

So this is the vacuum Rabi oscillation in the reducible representation, and the last term is the binomial distribution 
for N trials, with single-trial probability of success Z p . There are TV different frequencies and thus collapses and 
revivals will necessarily occur if 1 < TV < oo. For TV large enough the binomial distribution can be approximated 
by a Gaussian, and one can show that for small Z the parameter that controls the Rabi oscillation is effectively the 
product NZ (see below, Sec. VIII). The limit TV — > oo (with fixed Z) can be computed on the basis of the law of large 
numbers for the binomial distribution, 

r M 1 I 12 v sinVAV4+|ff|%t 

i.e. the frequency s/TV approaches the probability of success in a single trial of the Bernoulli process, s/TV — » Z p . Ijl4(l 
is essentially the standard Jaynes-Cummings prediction, but with a modified coupling. It is clear that the measurable 
coupling is not just g but rather its renormalized version g p ^ — g\fZ ' . Let us note that this is equivalent to bare 
charge renormalization: e D h = e§\fZ. Z is therefore an analogue of the renormalization constant Z% and Xfc = Z k /Z 
plays a role of a cut-off |27|. 

Both the cut-off and the renormalization constant occur here automatically. If we assume that for optical frequencies 
Z p = max k {Z k } = Z (i.e. Xp = 1) the agreement between the irreducible case and the TV — > oo limit of the reducible 
one is exact. The law of large numbers plays here a role of a correspondence principle with the standard formalism, 
a property not limited only to the Jaynes-Cummings example. 

With this background in mind one can easily generalize the discussion to thermal and coherent states, and mixed 
atomic initial condition [l5| . Let p + and p_ denote initial probabilities of finding the atom in, respectively, excited 
and ground states. Replacing vacuum by a thermal light with the distrubution 

n n 

PW = (1 + »)<"-») ■ (15) 



we find, p+(t) = w(t) + 1/2, 
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(16) 



s=0 x ' n=0 

The limit N — > oo 



^ V 1+nJ A 2 /4+\g ph \ 2 {n + l)xp 

is, up to Xpi known from irreducible representations. For a coherent state \z) and p+ = 1 we find 



„ M i y y l.9 P h| 2 (n + l) a sm 2 (t^Ay4+\ gph \2(n + l) S /(ZN)) ^/g^ .r^ ( N \ 

Z N AV4+W(n+lW(ZJV) »! V ?( P) 



The limiting form, for iV — > oo, is again familiar 



hm - l-f| n P(„ + lk sin 2 (VAV4 + | gph | 2 (n + l) Xp )|z phX ;|^ 



_ . ... n\ 

n=0 



(18) 



For the same reason as before we obtain the standard formula but with the cut-offs Xp = ^p/Z, x' P — y/Xp-i an d 
renormalized e p h = sq^/Z, z p ^ = zVZ. 

IV. DISSIPATIVE AND NONDISSIPATIVE DECOHERENCE 

An analysis of realistic experiments must take into account decoherence. There are two main sources of decoherence 
that were identified in the literature in the context of the experiment of Brune et al. [Tsj . 

The analysis of dissipation based on quantum trajectories approach [TtIEsI ] leads to the conclusion that the damping 
due to energy loss in the cavity should have the form p K , + (t) = e~ Kt p + (t), where p+(t) is the probability of finding 
the atom in the excited state in an ideal cavity, and 2n = 1/T cav . The factor 2 takes into account the fact that energy 
is not dissipated if the atom is in the excited state and there is no photon in the cavity. Obviously, for t — > oo the 
atom is with certainty found in its ground state. 

The second source of decoherence is nondissipative in nature and was discussed in . It originates from the fact 
that the data collected at time t should not be compared directly with pit) describing the state computed on the basis 
of first principles, but with the average 

/>oo 

PAt{t) = / dt'p At {t,t')p(t'), (19) 



Jo 

where pit') is the first-principles state and pAt(t, t') describes our lack of knowledge as to the exact duration of time 
evolution. The data from [lg involve a sample of 90 points selected from the time interval < t < 90 /its. Therefore 
the time-of-measurement uncertainty may be assumed to satisfy < At < 1 /is, which indeed turns out to reasonable 
model the data. However, it is not evident if this is really the true explanation of the discrepancy. The problem is 
that another value of At is also mentioned in ^(|, namely At = O.Oli. It would lead to a linear growth of At between 
0.01 /tzs and 0.9 /is, and then the agreement between theory and experiment is worse. 

Nevertheless, leaving aside this and similar subtleties, we may use the probability distribution introduced in 0] 

= At nt/At) ■ ( 20 ) 

In all the representations discussed in this paper we have arrived at atomic probabilities involving terms of the form 

p+(t) = A + B sin 2 fit. (21) 



The associated effective probabilities then read 
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FIG. 1: Standard theory of p+(t) for three different uncertainties of atomic time-of-flight measurements: At = 0.3 /is (dotted), 
At = 0.5 /is (full, suggested in At = 0.7 /is (dashed). Left plots: T cav = 220 /is, k = l/(2T cav ). Right plots: k = 0. 



POO 

PAt, K ,+ (t) = / dt' p At (t,t')e- Kt ' (A + B sin 2 SW) 
Jo 



A + -B11- 



2nAt 

1 + nAt 



is* /t 2QAt 
cos -t— arctan ■ 
VAt 



2QAt \ 
1 + rcAt/ 



(22) 



The overall damping factor (1 + nAt) t//At is the deformed exponential |29l| occurring in non-extensive thermodynamics 
|30| . and whose links to Gamma- function averages are well known |3l| . 



V. EXPERIMENT OF THE PARIS GROUP 



VACUUM RABI OSCILLATION 



Let us first concentrate on the vacuum Rabi oscillation observed in [18J. This part of the data is particularly 
intriguing and plays an important role for calibration of the experimental setup. Theoretical fits shown in 1|| were 
based on sinusoids exponentially damped by exp(— t/T) with T — 40 /is 32] . It is essential that the parameter T was 
much smaller from the reported value T cav = 220 /is of the cavity lifetime. The coupling constant employed in the fits 
was Sph/t = 47 kHz, and the cavity was filled with 0.8 K thermal light (average number of photons n = 0.05). Brunc 
et al. tried to explain the difference between T and T cav by means of dark counts and collisions with background gas. 

The role of dark counts and collisions was analyzed in detail in fl7i . but the conclusion was negative — the source 
of the discrepancy had to be different. The analysis presented in [TtJ revealed also another problem with the data: 
The cavity lifetime T cav = 220 /is should induce a shift of excited-state probability towards zero (as in Fig. 1, left 
part), but there is no trace of this phenomenon. 

Let us now turn to the solution of the discrepancy between T and T cav in terms of nondissipative decoherence, 
suggested in [l^. In Fig. I we show the prediction involving both kinds of decoherence, and based on irreducible 
representations of CCR. The initial probability is p+ = 0.99. For the left plots the damping factor is k = l/(2T cav ) = 
10 6 /440 Hz, and three different values of At are compared. The value At w 0.5 //s (suggested in would reasonably 
reproduce the data if one neglected the energy decay, a fact consistent with the observations from ■ I n the right 
plots k — 0. For times < t < 15 /is the data are then consistent with At = 0.7 /is. The next peak is well described 
by At = 0.3 fis, then again the minimum looks like At = 0.7 /is, and finally we can use At = 0.5 /is. The fits are quite 
sensitive to small variations of At, but generally the solution is acceptable if one could explain why k = is here 
meaningful. The problem is a serious one since the whole logic of the experiment is based on nonncgligible dissipation 
(cavity decay eliminates the maser efect). 



VI. EXPERIMENT VS. N < oo REPRESENTATIONS 



Let us now turn to the case of N < oo representations. The first question we have to clarify is what is the role of k 
for predictions based on reducible representations. In Fig. 2 (left) we compare the standard prediction from Fig. 1 for 
T cav = 220 /is and At = 0.5 /is (dotted) with an analogous result for the reducible N — 2000, Z = 0.1 representation 
(full). The two curves differ by less than experimental error bars, and it is clear that the finite- AT representations 
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p.(t) p,(t) 




FIG. 2: Comparison of the standard theory of p+(t) (dotted) with the reducible representation characterized by N — 2000, 
Z = 0.1 (full). The left plots, k = l/(2T cav ) = 10 6 /440 Hz, are shifted downwards with respect to the data. The right plots 
employ k = 0. All the curves correspond to At = 0.5 /xs. Predictions of the two theories differ by less than experimental error 
bars. 



suffer from the same problem as the irreducible ones: For n = l/(2T cav ) > the Lindblad-type plots are shifted 
downwards with respect to the data. This is not surprising, since for N oo the reducible representation should 
reconstruct predictions of the irreducible one. 

In Fig. 2 (right) we show the same situation as in Fig. 2 (left) but now with k = 0. The data are consistent with 
NZ = 200. It is very important to keep in mind that for N Z > 200 the agreement between the two theories will be 
even better. This is why this type of experiment will not be able produce an exact value of NZ, but only set a lower 
bound on the value of this parameter. Fig. 3 shows analogous plots for coherent states. Collecting all the avaliable 
data we can estimate a common lower bound following from various experimental situations — here in all the plots 
the lower bound NZ > 200 is enough to have predictions experimentally indistinguishable from the standard theory. 
Possibility of a test directly determining NZ is discussed in the next section. 



VII. CAN WE DIRECTLY MEASURE NZ? 



Reducible and irreducible representations are idistinguishable as long as the beats typical of finite N are masked 
by the decay caused by a nonzero At. In Fig. 4 we show the dynamics of p+(t) monitored with the resolution 
At = 0.005 /is. We assume that initially the atom is in the upper level and there are no photons (exact vacuum 
state at zero temperature). The plots reveal two important features of finite- N representations. First of all, even in 
exact vacuum we find beats analogous to what is known from irreducible-representation coherent states. Secondly, 
the first revival occurs after a time that depends effectively on the product NZ, and not separately on N and Z. To 
understand why this has to happen we replace the binomial distribution by its asymptotic form, valid for large N, 

(a-JVZ) 2 



N \ Ar / N \ nr e 2JVZ(1-Z) g IN Z 

Zi(\ - Z p f- S ~ Z s (l - Z) N - S » ; « - (23) 

3 J P P V 3 J y/2nNZ(l - Z) V2^NZ V ' 

The shape of the Gaussian is controlled mainly by the product NZ. The smaller Z the less important its exact 
value (the last approximate equality holds for small Z). Increasing N with Z kept constant we shift the first revival 
more to the right. In the limit N — *■ oo the first revival is shifted to infinity, and we recover the standard undamped 
oscillation. So, the absence of the revival in an experiment can only set a lower bound on NZ, and is not a proof that 
the physical representation is irreducible. 



VIII. FINAL REMARKS 



It would be interesting to analyze the other experiments involving finite- level atoms, especially those with masers 
and mazers [3(| l^l HH HH but a technical difficulty is that exact solutions are not there available at the moment. 
However, the experiments testing spectra of light |40l lilL l42l 14a. Ill ] in cavity QED are another realistic goal in this 
context. We have already computed the vacuum Rabi splitting, with the conclusion that for N — > oo we reconstruct 
the standard results, which is another example of the correspondence principle. The work on comparison of the theory 
with experiment is in progress, and we will present the results in a separate paper. 



FIG. 3: Comparison of the standard theory of p+(t) (dotted) with the reducible representation characterized by N — 2000, 
Z = 0.1 (full). Coherent states with n = 0.4 (left), n — 0.85 (right), and k — 0, At — 0.5 /is. Initially the atom is in the upper 
state (p+ = 1). 
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FIG. 4: Reducible representations for NZ = 200, n = 0, and At = 0.005 /is. Vacuum Rabi oscillation is monitored for a longer 
time. Left plot: N = 600, Z — 1/3. Right plot: N = 3000, Z = 1/15. For small Z and large N the revival occurs after a time 
that depends only on the single parameter NZ. 



The structure of vacuum collapses and revivals is like a fingerprint of the representation. The parameter NZ deter- 
mines the distance in time between the reviving peaks. For physical reasons Z must be a very small nonzero number, 
and thus N has to be very large, although finite. Confirmation that N < oo is physical would have fundamental 
consequences for renormalization theory, vacuum energy with all its implications, and studies of entanglement in 
cavity QED. The correspondence principle turns N < oo theories into generalizations of standard quantum optics. 
Our discussion explains why it is very unlikely that N < oo can be found inconsistent with experiment. And this is 
interesting in itself. 

We are indebted to M. Brune for the data, and Y.-T. Chough, D. R. Finkclstcin, S. Harochc, J. Naudts, G. Nogues, 
and W. Schleich for various comments. This work was done as a part of the Polish Ministry of Scientific Research 
and Information Technology (solicited) project PZB-MIN 008/P03/2003. 
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